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XMM- Newton observation of the Seyfert 1.8 ESO 113-G010: 
discovery of a highly redshifted iron line at 5.4 keV. 
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Abstract. We present a spectral analysis of the Seyfert 1.8 ESO 113-G010 observed with XMM-Newton for 4ks. 
The spectrum shows a soft excess below 0.7 keV and more interestingly a narrow emission Gaussian line at 5.4 
keV (in its rest-frame), most probably originating from a redshifted iron Ka line. No significant line at or above 
6.4 keV is found contrary to other objects showing redshifted lines, ruling out a strong blue-wing to the line 
profile. The line is detected at 99% confidence, from performing Monte Carlo simulations which fully account for 
the range of energies where a narrow iron line is likely to occur. The energy of the line could indicate emission 
from relativistic (0.17-0.23 c) ejected matter moving away from the observer, as proposed for Mrk 766 by Turner 
et al. (2004). Alternatively, the emission from a narrow annulus at the surface of the accretion disk is unlikely 
due to the very small inclination angle (i.e. less than 10°) required to explain the narrow, redshifted line in this 
intermediate Seyfert galaxy. However emission from a small, localized hot-spot on the disk, occurring within a 
fraction of a complete disk orbit, could also explain the redshifted line. This scenario would be directly testable 
in a longer observation, as one would see significant variations in the energy and intensity of the line within an 
orbital timescale. 
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1. Introduction 

In Active Galactic Nuclei (AGN), from Seyfert galaxies to 
quasars, the analysis of several X-ray features can help us 
to understand the central region of these powerful objects. 
Especially in the hard X-ray band, the FcKcv line complex 
observed in the 6-7 keV range is an important spectral di- 
agnostic tool to probe dense matter from the inner disk 
out to the Broad Line Region and the molecular torus 
(see review in Reynolds & Nowak 2003) . Recently, narrow 
spectral features in the 5-6 keV energy range were dis- 
covered with XMM-Newton and Chandra in a few AGN: 
NGC 3516 (Turner et al. 2002), ESO 198-G24 (Guainazzi 
2003, Bianchi et al. 2004), NGC 7314 (Yaqoob et al. 2003), 
Mrk 766 (Turner et al. 2004), and IC4329A (McKernan 
& Yaqoob 2004). Localized spots or narrow annuli which 
occur on the surface of an accretion disk following its il- 
lumination by flares has been proposed to explain these 
features (e.g., Nayakshin & Kazanas 2001, Turner et al. 
2002, Dovciak et al. 2004). An alternative scenario has 
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been proposed by Turner et al. (2004) for Mrk 766 for 
which the first evidence for a significant line energy shift 
has been observed over a few tens of ks. They proposed 
that this shift could be interpreted as deceleration of an 
ejected blob of gas traveling close to the escape velocity. 

Here, we present the first observation above 2 keV by 
XMM-Newton of the Seyfert type 1.8 galaxy ESO 113- 
G010 (z=0.0257, FWHM Ha = 2 000 km s^ 1 , Pietsch et 
al. 1998). The present XMM-Newton data shows that 
this object shows a soft excess and more interestingly an 
emission feature near 5.4 keV (in the Seyfert rest-frame). 
H =75kms _1 Mpc" 1 and qo=0.5 are assumed through- 
out. 

2. XMM-Newton observation 

ESO113-G010 was observed by XMM-Newton on May 

3, 2001 (OBSID:0103861601). The EPIC MOS cameras 
(Turner et al. 2001) were operating in partial window 
mode (exposure time about 6.9 ks), and the EPIC PN 
camera (Striider et al. 2001) was operating in full win- 
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dow mode (exposure time about 4.1 ks). The data were 
re-processed and cleaned using the XMM-Newton SAS 
version 5.4.1 (Science Analysis Software) package. Data 
are selected using event patterns 0-4 and 0-12 for pn and 
MOS, respectively, whilst only good X-ray events (with 
'FLAG=0') are included. The source spectra are extracted 
from a circular source region of about 32.6" radius, and 
the background spectra are extracted from a box (6' x 3') 
close to the source. The pn net exposure time is 4ks. The 
S/N of the MOS data is too low above 5 keV, therefore we 
use only the PN data which has the greatest sensitivity in 
the iron K-shell band. Unfortunately, for the present ex- 
posure time the S/N of the RGS data are by far too low 
for any spectral analysis. 

3. Spectral analysis 

In the following analysis, abundances are set to those 
of Anders & Grcvesse (1989), whilst we use the up- 
dated cross-sections for X-ray absorption by the inter- 
stellar medium from Wilms et al. (2000). The Galactic 
column density value is 2.74 x 10 20 cm~ 2 , obtained from 
the COLDENS program using the compilations of Dickey & 
Lockman (1990). We bin the background-subtracted spec- 
trum in order to have a minimum of 20 counts per bin. 
The errors quoted correspond to 90% confidence ranges 
for one interesting parameter (Ax 2 — 2.71). 

First we fit the data in the 1-4 keV energy range where 
the spectrum is relatively unaffected by the presence of 
a soft excess below about 0.7 keV, by a possible warm 
absorber, and by a Fe K emission line. Over this energy 
range, the data are well fitted by a single power law model 
with r=2.07±0.08 (x 2 =113.4/116). Figure 1 displays the 
spectrum extrapolated over the 0.3-8 keV broad band en- 
ergy. A strong positive residual is clearly seen below about 
0.7 keV due to the presence of a soft X-ray excess. In addi- 
tion, a clear positive deviation near 5.4 keV (in the Seyfert 
rest-frame) is observed. 



3.1. The soft X-ray spectrum 

The spectrum below 2 keV is only moderately fitted 
(% 2 /dof=268. 5/197) by an absorbed (by Galactic absorp- 
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Fig. 1. PN spectrum of ESO 113-G010 (z=0.0257, observer 
frame). A power law has been fitted to the 1-4 keV data and 
extrapolated to lower and higher energies. A soft X-ray excess 
is clearly seen extending to about 0.7 keV, as well as a positive 
deviation near 5.4 keV. 
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Fig. 2. A close-up of the data/model ratio in the 3-8 keV en- 
ergy range. The model is a Galactic absorbed power law. A 
clear positive residual is seen near 5.4 keV. 



body spectrum (diskpn in XSPEC, Gierlinski ct al. 1999) 
and a power law continuum, for the soft and hard band, 
respectively. We find a rather good fit with kT=121+ 1 [j eV 
and r=1.94±0.07 (x 2 /d.o.f=267.9/250). 
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(due mainly to the presence of a soft excess below 0.7 keV) . 
We find an upper limit for intrinsic absorption at the red- 
shift of the Seyfert of 0.8 x 10 20 cm~ 2 . Since this value is 
much weaker than the Galactic absorption, we fix the ab- 
sorption column density to that of the Galaxy (i.e., 2.74 x 
10 20 cm -2 ). We inferred a 0.1-2.4keV luminosity of about 
1.4 x 10 43 ergs _1 , which is compatible with the previous 
observations with ROSAT (Pietsch et al. 1998). The stan- 
dard explanation for the soft X-ray excess is that it results 
from thermal emission originating directly from the hot in- 
ner accretion disk (Malkan & Sargent 1982), and hence it 
is the high energy tail of the so-called Big Blue Bump. We 
fit the data with the combination of an accretion disk black 



Fitting the data above 1 keV with a (Galactic) ab- 
sorbed power law gives a good fit with r=1.98±0.06 
(x 2 /d.o.f=125. 1/128). The inferred unabsorbed 2-10keV 
flux is 2.4 x 10~ 12 ergcm~ 2 s _1 , corresponding to a lu- 
minosity of 3.1 x 10 42 ergs _1 . There is a clear positive 
residual near 5.4 keV as also shown in Figure 2. The fea- 
ture is well modeled by a Gaussian line (with cr=0.1keV) 
with A% 2 =8. 7 for 2 additional parameters. The line energy 
is 5.38±0.11keV with EW=265± 1 4 JeV. Figure 3 shows 
the confidence contour plot for the rest-energy and inten- 
sity of the line corresponding to the 68%, 90% and 95% 
confidence levels. One should notice that a Cr Ka (at 
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5.4 keV) is present in the internal EPIC-pn background 
(see Figure 2 in Freyberg et al. 2002) but its strength 
is negligible and therefore can be excluded as a possibil- 
ity to explain the 5.4 keV line. Formally, we do not find 
any statistical requirement for a narrow line at 6.4 keV 
(Ax 2 =0.8 for 1 additional parameter). Fixing the line en- 
ergy at 6.4 keV and the width at 10 eV (i.e an intrinsically 
narrow line), we find an upper limit for the EW of 230 eV. 
Splitting the observation into two equal parts (of about 
2ks each), we do not find any significant line parameter 
variation during this very short timescale. 



3.3. Estimating the detection probability of the line 

In estimating the significance of the 5.4 keV line, straight- 
forward application of the standard two parameter F-test 
does not take into account the range of energies where 
lines can occur in the spectrum, nor does it easily account 
for the number of bins in which one is searching for lines. 
The probability P of finding a feature at one particu- 
lar given energy can be estimated from the F-test, upon 
the addition of one extra degree of freedom (the line nor- 
malisation), which in ESO 113-G010 is 99.8%. However in 
searching for rcdshiftcd (or blueshifted) iron lines, we are 
fitting over many bins or resolution elements, where nar- 
row features may occur by chance (i.e. statistical noise). 
Thus the probability of detecting a line at any energy 
in the iron K-shell band can be estimated by Pjv = Pq 1 , 
where N is the number of bins searched. In ESO 113-G010, 
we define the iron K band between 4-7 keV (i.e. where one 
may realistically expect to find iron lines), which consists 
of 18 spectral bins (at 20 counts per bin). Thus the prob- 
ability of detection of the 5.4 keV feature by this estimate 
is 96.5%. Note this is probably an under-estimate of the 
detection probability, as some of the actual spectral bins 
are smaller than the instrument resolution and therefore 
the number of bins (or trials) where a line may be located 
has probably been over-estimated. 



We also carried out a more rigorous test of the signif- 
icance of the line using Monte Carlo simulations. For our 
null hypothesis, we assumed that the spectrum is simply 
an absorbed power-law continuum, with the same param- 
eters as the absorbed power-law model fitted to the real 
data. We used the xspec fakeit command to create 1000 
fake EPIC-pn spectra corresponding to this model, with 
photon statistics expected from the 4 ks exposure, and 
grouped each spectrum to a maximum of 20 counts per 
bin. Following the procedure used to test the real data for 
the presence of a narrow line, we fitted each fake spectrum 
with an absorbed power-law (absorption fixed at Galactic, 
but power-law photon index and normalisation left free to 
vary), to obtain a x 2 value. We then added a narrow line 
(<r = 0.1 keV) to the fit, restricting the line energy to be 
between 4-7 keV. Furthermore we stepped the line over 
the 4-7 keV energy range in increments of 0.1 keV, whilst 
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Fig. 3. Contour plot showing the 68%, 90% and 95% (from the 
inner to outer curves) confidence level for values of line rest- 
energy (keV) and line flux for the line detected near 5.4 keV. 
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Fig. 4. The measured change in x 2 against cumulative prob- 
ability for the addition of a single Gaussian line in 1000 ran- 
domly generated simulated spectra. The null hypothesis for the 
simulated spectra is a power-law continuum with Galactic ab- 
sorption. For ESO 113-G010, the measured Ax 2 = 8.7 for the 
addition of a Gaussian line at 5.4 keV is only exceeded in 1.4% 
of the simulated spectra. Therefore the 5.4 keV line is detected 
at ~ 99% confidence. 



fitting separately each time to ensure the lowest \ 2 value 
was found. We then recorded the minimum \ 2 obtained 
from these multiple line fits for each fake spectrum, and 
compared with the corresponding x 2 01 the null hypothe- 
sis fits, to obtain 1000 simulated values of the Ax 2 , which 
we used to construct a cumulative frequency distribution 
of the Ax 2 expected for a blind line search in the 4-7 keV 
range, assuming the null hypothesis of a simple power- 
law with no line is correct. The resulting distribution of 
X 2 against detection probability is shown in Figure 4. We 
therefore find only 1.4% of fake power-law spectra fitted 
with a line show a larger Ax 2 than observed in the real 
data, implying that the line detection is significant at ap- 
proximately 99% confidence. 
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Fig. 5. Line profiles obtained by a relativistic line model for 
two inclination angles (9). Top: #=5°. Bottom: 9=20°. 

3.4. A relativistic disk origin for the line? 

The 5.4 keV line could be a relativistic line for which we 
only see one horn due to the moderate spectral resolu- 
tion. Therefore, we fit the data with a diskline profile 
(Fabian et al. 1989). We assume a narrow annulus of 1 
R g , in order to obtain a narrow profile. Initially we fix 
the inclination angle at 9 — 45°, and assume the line is 
emitted at 6.4 keV (i.e near-neutral iron). However this 
model fails to re-produce the narrow emission line at 5.4 
keV. Upon varying the disk inclination angle a good fit to 
the line is obtained, with x 2 / d -°- f = n 8-4/127, requiring 
a very small inclination angle of 9 ~ 5°, with a formal 
upper- limit at 90% confidence of 9 < 10°. We constrain 
the inner radius of the annulus at Rj n =9.5lg;g Rg, whilst 
the line EW is 303±166cV. Thus for a very low (face-on) 
disk inclination the line is gravitationally redshifted whilst 
there is very little transverse velocity shift, hence the red 
and blue wings of the line are not resolved. A narrow red- 
shifted feature is then observed, consistent with the line 
profile in Figure 2. At higher disk inclination angles, the 
transverse velocity shift is much more prominent leading 
to a strong blue wing, which is not observed in the data. 
This is illustrated in Figure 5 where the theoretical line 
profiles at two angle inclinations are compared, i.e., 9 =5° 
and (9=20°. 

4. Discussion 

It has been proposed by Skibo (1997) that energetic pro- 
tons can be responsible of the destruction (spallation) of 
Fe into Cr and other lower Z metals on the surface of the 
accretion disk, which enhances the line emission expected 



from elements of low abundance. The strongest line apart 
the FeKa line expected is the CrKa line at 5.4 keV, and 
the second one being the MnKa at 5.9 keV (see Figure 
4 in Skibo 1997). The energy of the feature we find here 
is compatible with the energy of the CrKa line. However 
the ratio of the line at 5.4 keV (if CrKa) and the upper 
limit of a narrow Gaussian at 6.4 keV (due to iron) is not 
consistent with the line ratio predicted by this model, i.e. 
the line flux of the CrKa should represent 33% of the 
FeKa. This is not compatible with EW=265t \f 5 eV, and 
EW<230eV found respectively at 5.4 keV and 6.4 keV. 
Therefore iron spallation seems to be an unlikely expla- 
nation in this also found for other similar line 
observations (e.g., NGC3516 and Mkn766). 

If instead the line originates from an annulus on the 
disk, then the small inclination angle derived above does 
not appear to be consistent with the classification of this 
object as a type 1.8 Seyfert, as one might expect a larger 
inclination angle of about 45-60°. However, Nayakshin & 
Kazanas 2001, Turner et al. (2002) for NGC3516, and 
Dovciak et al. (2004) have proposed that these features 
could be due to localized hotspots on the accretion disk 
surface following its illumination by flares. In this case, the 
range of transverse velocity shifts would be much smaller 
than for an entire disk orbit, resulting in a relatively nar- 
row emission line feature, even at lower inclination angles. 
Narrow red-shifted features could be produced with an 
inclination angle of up to about 30 degrees in this model 
(see Fig. 1 in Nayakshin & Kasanas 2001, and also Fig. 3 
in Dovciak et al. 2004), while at higher inclination an- 
gles the line becomes blue-shifted above 6.4 keV, which 
is ruled out from the data. The narrow width of the ob- 
served emission line would imply that the emitting region 
must be small and that it is observed for only a fraction 
of the entire orbit. Indeed the orbital timescale of matter 
co-rotating along a circular trajectory assuming r=10i? g is 
between 2.7-27 hr for a black mass of 10 7 M Q and 10 s M , 
respectively, probably larger than the exposure time of 
the present observation. A longer exposure is therefore re- 
quired to study possible variations in the line profile as a 
function of phase of the spot. 

The redshifted line emission could be also due to the 
signature of an inflow or outflow of matter from mate- 
rial close to the nucleus. Moving material that is not part 
of the disk structure would be sensitive to the effects of 
strong gravity and may achieve velocities that are a signif- 
icant fraction of the speed of light, causing significant dis- 
placement of emission lines from their rest energy. In the 
XMM-Newton spectrum of Mrk766, Turner et al. (2004) 
found an energy shift from 5.60 to 5.75 keV in a line oc- 
curring over a few tens of ks. They proposed that such 
rapid energy shift could be interpreted as deceleration of 
an outflowing ejected blob of gas traveling to the escape 
velocity, where the material is viewed on the reverse side 
of the black hole (i.e. moving away from the observer). 
The energy of the line in ESO 113-G010 could indicate a 
(redshifted) fast moving medium at 0. 16-0.23 c, assuming 



Porquet et al.: XMM-Newton observation of ESO 113-G010 



5 



a neutral iron (6.4 keV in its rest-frame) or a H-like iron 
(7keV in its rest-frame) respectively. 

Up to now mainly red-shifted narrow lines (except for 
the NLS1 RX J0136. 9-3510, Ghosh et al. 2004) have been 
observed in Seyfcrt 1 galaxies. For the hotspot model, 
this could be explained by the fact that red-shifted lines 
are predicted for very small disk inclination (up to about 
30 degree), which is compatible with the Unified Scheme 
(Antonucci 1993). For the blob model, this could be ex- 
plained also by a geometry effect. The outflow could be 
optically thick on the near side of the flow, so one only sees 
emission from the reverse side. However given the statisti- 
cal quality of this dataset at energies above the iron K line 
(> 7keV), it is not possible to exclude a blue-shifted com- 
ponent to the line from an outflow, which has the same 
strength as the redshifted 5.4 keV line. For instance the 
90% confidence upper limit on the equivalent width of a 
line detected between 7-8 keV is < 340 eV, i.e. consistent 
with the line observed at 5.4 keV. 

The present observation of ESO113-G010 is far too 
short (4 ks) to study any shift of the line energy and mea- 
sure any outflow or inflow of the material. Indeed within 
the current snapshot observation, three scenarios may be 
possible; (i) the line emitting matter may be accelerating 
due to gravity onto the accretion disk or black hole; (ii) the 
material may be ejected from the accretion disk and decel- 
erate under gravity away from the observer or (iii) the line 
may originate from a rapidly variable hot-spot on the disk 
surface. In case of a material accelerating or decelerating 
under gravity, we can use Av/ At = GM/R 2 (where At is 
the time under which the deceleration is observed and R is 
the radius at which the material is situated at that time) 
to estimate the observed time of velocity change. For ex- 
ample, assuming a change of the energy line from 5.4 keV 
to 5.6 keV, and R~ 3.5 x 10 14 cm (as for Mrk 766, Turner 
et al. 2004), we inferred At <~ 90 ks and At ~ 9ks, assum- 
ing a black hole mass of 10 7 M© and 10 8 M©, respectively. 
In the case of the hotspot scenario, both line energy and 
line intensity are predicted to vary with time on the or- 
bital timescale (from 10 ks up to 100 ks if the black hole 
mass is 10 7 M© and 10 8 M Q , respectively), as shown in 
Fig. 1 in Nayakshin & Kazanas (2001) (see also Figs. 3 
and 4 in Dovciak et al. 2004). Therefore a much longer 
XMM-Newton observation, of about 100 ks in length, is 
required to be able to differentiate between these models. 

In this paper we have performed Monte Carlo sim- 
ulations to calculate a realistic detection probability for 
the redshifted line at 5.4 keV. The Monte Carlo method 
takes into account that an iron line-like residual may be 
detected at any energy within the iron line band (defined 
as 4-7 keV in this instance); in contrast the conventional 
F-test is likely to over-estimate the detection significance, 
as this does not account for the range of energies or bins 
where a redshifted line may be observed. Applying the 
Monte-Carlo method to ESO113-G010 results in a de- 
tection significance of 99%, which is still just under the 
3cr confidence level. Clearly, a longer observation will be 



needed to increase the detection significance of the fea- 
ture to > 99% confidence. This is also the first time a 
Monte-Carlo probability estimate has been reported for 
one of the several previously published narrow redshifted 
line features in AGN, e.g., NGC 3516 (Turner et al. 2002), 
ESO 198-G24 (Guainazzi 2003, Bianchi et al. 2004), NGC 
7314 (Yaqoob et al. 2003), Mrk 766 (Turner et al. 2004), 
and IC 4329 A (McKernan & Yaqoob 2004). It would there- 
fore be desirable to re-evaluate earlier claims on redshifted 
line detections using more robust Monte-Carlo techniques, 
as well as analysing a wider sample of Seyfcrt galaxies to 
determine the frequency of the relativistically shifted lines 
in AGN. 
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